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Abstract: At the end of inflation, the coherent oscillations of the inflaton field may 
resonantly amplify the long wavelength modes of both bosons and fermions coupled 
to it. We study the resonant production of both kinds of particles during preheating 
in a model of hybrid inflation. The coherent time evolution of the inflaton and the 
Higgs fields after inflation induce a very different production of fermions depending 
on whether they are coupled to the Higgs or to the inflaton. For reasonable values of 
the model parameters, the fermion production through parametric resonance can be 
very efficient. We study the relative growth of the fermion and boson energy densities 
during preheating in hybrid models. During the initial stage of preheating, fermion 
production dominates the relative energy density, while the exponential growth of 
bosonic modes soon takes over. 
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1. Introduction 

Cosmological inflation is an extremely efficient mechanism in diluting any particle 
species or fluctuations. At the end of inflation, the universe is empty and extremely 
cold, dominated by the homogeneous coherent mode of the inflaton. Its potential 
energy density is converted into particles, as the inflaton field oscillates around the 
minimum of its potential. These particles are initially very far from equilibrium, but 
they interact among themselves (or decay) and thermal equilibrium is achieved at a 
very large temperature. From there on the universe expanded isoentropically, cooling 
down as it expanded, in the way described by the standard hot Big Bang model. 
Recent developments in the theory of reheating suggest that the decay of the inflaton 
energy could have been explosive [|]], due to the coherent oscillations of the inflaton 
zero mode, which induce its stimulated decay. The result is a resonant production of 
particles in just a few inflaton oscillations, in a process known as preheating M. The 
number of particles produced in this way is exponentially large, which may account 
for the extraordinarily large entropy, of order 10 89 , in our observable patch of the 
universe today. 

Preheating is not generic, it may occur in different models of inflation p], |3|. Q 
but only under special circumstances. Preheating strongly depends on the inflaton 
couplings to other fields, as well as on the amplitude and frequency of the inflaton 
oscillations. This explosive production occurs through parametric resonance of the 
long wavelength modes of any field coupled to the inflaton, either bosonic |I| or 
fermionic || . Due to the parametric nature of the resonance the particle production 
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occurs for well defined frequency bands. The difference between bosonic or fermionic 
preheating is spectral: bosons can have large occupation numbers for a single mode, 
while fermions saturate due of the Pauli exclusion principle, and therefore a larger 
fraction of their energy can be transferred to higher momentum modes. Also, the 
rate of growth of these two components is very different: while the boson energy 
density grows exponentially []IJ, that of fermions cannot 0. As a consequence, very 
soon bosons dominate the decay of the inflaton energy if both kinds of particles are 
resonantly produced. After a few inflaton oscillations, the energy density in fermions 
and bosons may be large enough to backreact on the inflaton oscillations and hence 
this will eventually stop their resonant production [§, 

After preheating in chaotic inflation, the main part of the energy density of 
the universe is in the coherent modes of bosons: the inflaton, as well as any other 
parametrically amplified bosonic fields with large occupation numbers, while a sub- 
dominant part may be in the parametrically produced fermions. This state is very far 
from thermal equilibrium (characterized by a decoherent ensemble with small occu- 
pation numbers for all momentum modes) and, in fact, such a state could be used to 
generate the required baryon asymmetry during preheating, either via GUT baryoge- 
nesis || , EW baryogenesis or leptogenesis [pj . At the end, the final thermalization 
of the universe occurs through the rescattering of all the particles, which breaks the 
coherence of the bosonic fields. See, for instance, Ref. for thermalization via 
fermionic modes. 

In this paper we study the resonant production of both fermions and bosons in a 
hybrid model of inflation [|11|]. In these models it is possible to choose, without fine- 
tuning, the masses and couplings of the fields in such a way that the rate of expansion 



is negligible compared to the masses involved |C|. In that case, preheating occurs 
in a few inflaton oscillations after the end of inflation, before the scale factor has 
grown by a single e-fold, and therefore we can ignore the universe expansion. This 
reduces the problem to resonant production of fermions in Minkowski space, for 
which there is a complete analytic treatment ||. We find that in the case of hybrid 
inflation, contrary to what happens in chaotic inflation, it is possible to produce, 
for certain couplings of the fermions to the inflaton field, a larger contribution of 
fermionic modes than bosonic ones to the energy density of the universe at the end 
of preheating. 

We describe in Section 2 the hybrid model under consideration, with a Higgs- 
type field coupled to the inflaton, and discuss the solution of the classical equations 
of motion. In Section 3 we study the pair production of fermions in this model, for 
a non-zero coupling of fermions to both the inflaton and the Higgs field. We analyze 
the growth of the fermion energy density for different values of the couplings, as 
well as for different fermion masses. In Section 4 we study the boson production for 
both couplings to inflaton and Higgs, together with the growth of the bosonic energy 
density. In Section 5 we draw our conclusions. 
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2. Hybrid model 



We will consider a model of inflation with two fields, one that slow-rolls down a 
potential, driving inflation, and another one with a symmetry breaking potential 
that triggers the end of inflation, which we will call the Higgs. Such a model is 



called hybrid inflation and was proposed by Linde long ago [fUJl . There are many 
particle physics realizations of this general class of models fl3| , [14|, [15| . One particular 
potential for hybrid inflation is given by the tree level potential 1 

V(cf>, a) = ^- \m 2 <? + \xtf + \g 2 <?o 2 + imV . (2.1) 

During hybrid inflation, the inflaton field a evolves along the Higgs valley at = 0. 
As soon as the Higgs acquires a negative mass term, it triggers the end of inflation. 
That is, in less than one e-fold the two fields start oscillating around their common 
absolute minimum ((f) — v, a — 0), where v = Mj y/X is the Higgs vacuum expectation 
value (vev). The end of inflation occurs when the mass of the Higgs vanishes, i.e. at 
a = a c = M/g. 

By defining y = ojo c and x = <f)/v, we can write the potential ( |2.1|) as 

V(x, y) = V [(l - x 2 ) 2 + 2x 2 y 2 + 2 iy 2 } , 

where V = M 4 /AX, and 7 = Xm 2 /g 2 M 2 is a constant that is constrained to be 
small by the amplitude of the microwave background (CMB) anisotropies. It does 
not play any significant role after inflation and we will neglect it here. Also, the rate 
of expansion can be made relatively small in hybrid inflation (as long asM< gMp) 
and still satisfy the CMB constraints, see Ref . ]12[ . Therefore, we will ignore the rate 



of expansion of the universe here. We will also concentrate in the regime A ^> g 2 , for 
which the Higgs evolves along the minimum of the potential, following the inflaton 
oscillations. 

We can then write the evolution equations after inflation, redefining the time 
unit as x e Mt (with M = gM/VX), as: 

y" + (l-y 2 )y = 0, (2.2) 
l-y 2 = x 2 . (2.3) 

For initial conditions y(0) = yo and y'{0) = 0, we find the solution 

y = y cd(u\m) , 
u = r^a 2 /2 , m = ai/a 2 , (2.4) 

Oi = 2/o , 0,2 = 2 - yl , 



1 See Ref. [19 for a discussion of loop-corrections to this potential during inflation. 
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where cd(u|ra) is the Jacobi elliptic function. This solution is periodic in the r 
variable, with period 

T = 4* —K(m), (2.5) 

V a 2 

with K the elliptic integral. We have plotted the evolution of both the normalized 
Higgs and the inflaton after inflation in Fig. |I[ Throughout the paper we will take 
yo = 0.9999. The corresponding period of the inflaton oscillations is T ~ 29.97. 
Although this period increases with yo approaching one, the mean energy density 
transferred to fermions coupled to the inflaton is nearly unchanged (provided that 
yo is close to one). If the fermions are coupled to the Higgs, the results are more 
sensitive to the value of yo for strong coupling, as we will discuss later. 




0.5 1 1.5 2 2.5 

t/T 

Figure 1: The time evolution of the inflaton and Higgs fields after inflation, for initial 
condition yo = 0.9999. 



3. Pair production of fermions 

The oscillations of the inflaton and Higgs fields at the end of inflation trigger the 
explosive production of particles that couple to either or both of these fields. The 
production of bosons in hybrid inflation has been studied in Ref . p, [H] . Here we will 
analyze the fermion production in this model and in the next Section will elaborate 
further on the bosonic case in order to compare the two results. 

We will consider first the coupling of fermions to the inflaton a, with coupling 
hiaipip, and a possible mass term m^ipip. Then we will consider the coupling of 
fermions to the Higgs, with coupling h2<pif)ip- When the symmetry gets broken this 
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last coupling will give the fermions a mass through the non vanishing vev of the 
Higgs. 



3.1 Coupling to the inflaton 

Let us consider here a fermionic field ip satisfying the Dirac equation 

- ha(t) - m f ) if> = 0. (3.1) 

The solutions are more easily obtained using an auxiliary field X(x,t), such that 
ip = (ij^Vfi + hia(t) + m^) X. Decomposing it as exp(iK ■ x)X K (t)R r , with R r 
eigenvectors of 7 with eigenvalue +1, we can write the equation of motion for fermion 
modes X k as 

X£ + n 2 k X k -iy/qiy , X k = 0, (3.2) 

where 

tt 2 k (T) = k 2 + (^- iy (r)+m^) 2 , (3.3) 

fc^^S (3-4) 
ifif = m^/M (3.5) 

and we have rescaled k = K/M. We will here display the results for qi between 1 
and 10 6 . Notice that if we were to take e.g. A = 1 and g = 0.01 this will lead to 
qi = 10 8 h 2 , and hence the values of q\ considered would correspond to h\ between 
10" 4 and 0.1. 

We take the initial conditions corresponding to positive frequency plane waves 
at t < 0: 

X k (0) = [2tt k (tt k + y^yo + m^)}- 1 ' 2 , (3.6) 
X' k (0) = -iQ k X k (0). (3.7) 

The quantity of interest to us is the fraction of the total energy, p to tai — Kb 
which is transferred into fermions, 

— = h\^- \dkk 2 Q k (r)n k (r) , (3.8) 

with the number of fermion-pair degrees of freedom iV = 1 or 2 for Major ana or 
Dirac fields, respectively. The occupation number (for fermion pairs), n k , can be 
calculated as JTS[ 



Z il k Z\L k 
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Figure 2: The time evolution of the fermion occupation number n k (continuous line) and 
the averaged particle number n k (dashed line), for values of k 2 = 0.51 and 1.4, correspond- 
ing to the first peaks of F k for q\ = 1 and 100 respectively. 



Note that n k (0) = 0, thanks to the initial conditions ( |3.6|) and (|3.7| ). Also, it is easy 
to see that the occupation number of fermion pairs is always smaller than one, 2 as 
expected. 

Particle production through parametric resonance occurs at those moments when 
the adiabaticity condition dQ k /dr < VL 2 k is violated, see Ref. 0. For fermion pro- 
duction this will occur in general whenever the effective fermion mass rh^ + y/qiy 
approaches zero. This is illustrated in Fig. |2|, where the continuous lines show the 
evolution of the exact occupation number for massless fermions and two values of the 
qi parameter. The occupation numbers n k (solid lines) have jumps every quarter and 
three quarters of the irmaton period, corresponding to the times when the irmaton 
value vanishes, y = (see Fig. H), and therefore when the effective fermion mass also 
vanishes. 

To evaluate Eq. ( |3.8|) using n k {r) from Eq. (|3.9|) and the numerical solutions of 
( |3.2|) is however quite demanding, so that it is convenient to use the approximate 
analytical method developed in [18, [5]. This method exploits the periodicity of y 
and allows to obtain a smooth function 



n k (r) 



F fc sin (u k r) 



(3.10) 



which coincides with n k (r) for r = nT. Hence, n k (r) gives an approximate expression 
for n k (r) without its fine (spiky) details. The approximate solution (dashed lines in 

< 1, for fermions, in terms of the 



2 This is related to the fact that n k — |/3fe| 2 = 1 — \a k \ 2 
Bogoliubov transformation coefficients {a k ,f3k), see Ref. [[L8[. 
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Fig. HD follows the overall oscillations of the occupation number, matching the exact 
results at every inflaton period (and also at every half period in this case). 




Figure 3: The instability chart for fermions coupled to the inflaton. Contours correspond 
to equipotential values of Fk=0.1, 0.5 and 0.9 in the plane (qi,k 2 ) from lighter to darker. 

The advantage of using fifc(r) is that it has a simple temporal behaviour, while 
its k dependence can be obtained from the knowledge of the functions Fj. and u^, 
where 

F k 

cos (z/ fc T) 

and Xjp satisfies the same equation ( |3.2j ) with the initial condition X^^O) = 1, 
X^ (0) = 0. Therefore, to obtain n^r) we only need to solve Eq. ( |3.2|) during one 
inflaton period. The two signs in Eq. (|3.12|) correspond to two possible functions Uk, 
which oscillate with the same amplitude but different frequency, and both match 
at every period of the inflaton oscillation. The best approximation to n^ir) is given 
by the one with smaller frequency for the k 2 values where F& is maximum (i.e. 
for the momenta contributing significantly to p F ), which in this case corresponds to 
taking the minus sign in Eq. ( |3.12| ). This is the function plotted in Fig. Notice 
that the maximum value of v k is n/T ~ 0.105, as can be deduced from Eq. ( |3.12j ). 

We show in Fig. |3| the instability chart for massless fermions coupled to the 
inflaton. It displays the contours in the (qi,k 2 ) plane of equal F%. values. Fermions 



fc 2 (ImXW(T)) 2 
Q 2 (T) sin> fc T) 



±ReX«(T). 



(3.11) 
(3.12) 
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are mainly produced with momenta in the darker regions, corresponding to maxima 
of Fk- We see that the bands get narrower with increasing k 2 for a given q\ value, 
and after several bands they shrink to a negligible width. 

The upper panel in Fig. f| shows the behaviour of Fk for different q\ values 
extending up to q\ = 10 6 . These correspond to cuts in the instability chart (Fig. |3|) 
at a fixed value of q\. The maximum momentum for which the bands are sizeable 
grows as /c max ~ Qi^ for fermions coupled to the inflaton. 3 The lower panel of Fig. § 

3 Note that the same behaviour was found in Ref. M for chaotic inflation. However, a different 
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shows Vki as a function of k 2 , for the same values of q%. Note that the maxima of F k 
correspond to local minima of v k . 




2 4 6 8 10 

t/T 

Figure 5: The fraction of the total energy transferred to fermions coupled to the inflaton, 
as a function of time, for q\ = 1, 10 2 , 10 4 , 10 6 . 

The fraction of the total energy transferred to the fermions then finally results 
^l = hl^- [ dkk 2 n k (r)F k sm 2 (u k r). (3.13) 

Ptotal 7T Si J 

Its time evolution is shown in Fig. || for different values of the coupling parameter 
q\. We see that after the first oscillation of the inflaton this fraction already reaches 
its asymptotic value and then fluctuates around it. The asymptotic value p F /ptotai 
scales as ql^h 2 . Moreover, for large values of q± a significant fraction of the inflaton 
energy can be transferred into fermions (as long as hi is not too small). We have 
checked that the final density is insensitive to the initial value of the inflaton field 
y as long as \yo — 1| < 1CT 2 . This can be simply understood since, for = 0, the 
production of fermions takes place as the inflaton crosses through y = 0. 

If fermions are massive, the inflaton energy is transferred to the fermions with 
nearly the same efficiency as for massless fermions up to a maximum cut off value 
of the fermion's mass, above which it drastically drops. Fig. |6| shows this behaviour 
for different values of the q\ parameter. The cut off value of the mass goes like 
^i/),cutoff ~ Si^ 2 - 4 The reason is simple, particle production occurs whenever the ef- 
fective fermion mass vanishes, m eff = + \Zqiy(r) = 0. For small values of q±, even 

behaviour, fc max ~ q 1 ^ 3 , was found in Ref. Q and was attributed, in their case, to the redshift of 
the modes during the particle production process. 

4 The same behaviour was found for chaotic inflation in Ref. ||7| . 
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Figure 6: The fraction of the average energy density in fermions, h l 2 (p F ) / Ptot&\ as a 
function of the fermion mass (averaged between the eighth and tenth inflaton periods). 

a small (positive) rh^ will prevent m e g from vanishing as the inflaton y{r) oscillates. 
However, as we increase qi, larger bare masses are still allowed for particle produc- 
tion. The largest (negative) value of y(r) ~ —1 gives the cutoff mass. Recalling that 
the actual mass of the fermion is = (g/VTlm^M, and that we are working on 
the regime M < gMp, we see that this still allows for quite large fermion masses to 
be produced. 

3.2 Coupling to the Higgs 

We will consider now the coupling of fermions to the Higgs, and study their pair 
production. The associated Mathieu equation is analogous to fl3.2|) , but with x(r) = 
+yl — y(T~) 2 in the place of y(r), i.e. 

X' k ' + VL\X k - i^q- 2 x'X k = , (3.14) 

where Vt\ = k 2 + q 2 x 2 , x' = —y'y/x, and q 2 = h\/ g 2 . 

The initial conditions are given by (|3.6|) with y Q — > xq and we take the bare 
mass fh^ = 0, since we are assuming that the fermion acquires a mass through the 
Higgs mechanism. 

The occupation number is given by 

n k (T) = l-^Im[X k Xl*]-y^. (3.15) 
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(x - T/4) / T 



Figure 7: The time evolution of the fermion occupation number (continuous line) 
and the approximating function (dashed line), for values of k 2 = 0.235 and 0.371, 
corresponding to peaks of for c/2 = 1 and 100 respectively. 




2 4 6 8 10 



Figure 8: The instability chart for fermions coupled to the Higgs. Contours correspond 
to equipotential values of Ffc=0.1, 0.5 and 0.9 in the plane (q2,k 2 ). 

Particle production again occurs periodically. In this case, the jumps in the occu- 
pation number take place for every period and half period of the irmaton, i.e. at 
t = nT/2. These correspond to the times when the Higgs approaches zero, and thus 
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when the effective mass of the fermion is at a minimum. An approximate expression 
rtfc(r) can be constructed as in the previous case of fermions coupled to the inflaton. 
However, a function that matches n k {r) at the times when it has the spikes would 
not generally be a good approximation to the n^r). Thus we have performed a shift 
of a quarter of period in the initial time, so that n^fY) matches the value of n^r) 
at r = (2n + l)T/4, i.e. where n^r) has the plateau, and this gives a very good 
approximation as it is shown in Fig. [7]. Starting with n^{r) = at r = T/4 does not 
affect the results after a few inflaton oscillations, since the density rapidly saturates 
to its asymptotic value. 

The expressions for F k and v k are given by Eqs. (|3.11|) and ( 3.12|) , respectively. 
In this case the function with smaller frequency in the peaks of F k corresponds to 
the solution of Eq. (|3.12|) with the plus sign choice. We show in Fig. § the instability 
chart for fermions coupled to the Higgs. The resonance bands are narrower than 
those corresponding to the coupling to the inflaton. 

The upper panel in Fig. |S] shows Fk as a function of k 2 for several values of q 2 , 
between 1 and 10 4 , while the lower panel shows the corresponding v k . Note how 
quickly the bands become very narrow as k 2 is increased. As a consequence, the 
fermion production through the coupling to the Higgs is in general less efficient than 
through the coupling to the inflaton in hybrid inflation models. 

The ratio of Higgs-coupled fermion energy density to total energy is given by 



pjj) _ g 2 h\ 2N 

Ptotal A 7T 2 q 2 



dkk 2 n k (r)n k (T) . (3.16) 



and it is shown in Fig. [10]. If we take as a crude fit of the numerical results that 
(\/g 2 h 2 ,)p F /p total ~ KT 2 /g 2 , we see that p F /p to tai ~ 10~ 2 ^ 4 / A, and since we are 
working in the regime A ^> g 2 , we see that the fraction of energy transferred to 
fermions here remains small (but not necessarily negligible). 

Since the fermionic production takes place when x ~ 0, it is necessary to take the 
initial condition y sufficiently close to one so that small Higgs values occur. Taking, 
as we did, yo = 0.9999, ensures that at least x ~ 10 -2 is reached, and hence the 
parametric resonance is not inhibited by a mass gap for the values of q 2 considered. 

Notice that when the symmetry is finally broken, the fermions get a mass = 
h 2 M/ V% which also in this case can be quite large. 



4. Boson production 

The parametric resonant production of Higgs field particles has been shown to be 
quite inefficient in hybrid inflation models 0. However, other scalar particles coupled 
to the inflaton or the Higgs can can be abundantly produced. Here we will compute 
the parametric production of these bosons in order to compare it with the fermionic 
production. 
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Figure 9: The fermion (envelope) spectrum Fj~ (upper panel) and Uk (lower panel) as a 
function of k 2 , for fermions coupled to the Higgs, with qi = 1, 10 2 , 10 4 . Note how narrow 
the peaks become for increasing values of qi- 

We consider a scalar field x with mass m x and coupling to and a given by 

v(x) - \™\x* + ^ixV + (4.1) 

We can write the equation of motion for the scalar field modes X k as 

X£ + ^(r)X fc (r)=0, (4.2) 

where 

Ql(r) = k 2 + m 2 x + q [ y\r) + q' 2 (1 - y\r)) , (4.3) 
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/ X 9l ,92 m X (A A\ 

Q^ = — t , Oo = — — , m v = . 4.4 

1 g 4 g 2 M K ' 

We chose initial conditions of positive frequency plane waves at r < 

X k (0) = [2ft fc r 1/2 , (4.5) 

X' k (0) = -m k X k (0) . (4.6) 

The boson occupation number can then be calculated as 

•*M = ^i*ir+£w-5- < 4 - 7 ' 



Note that n k (0) = 0, thanks to Eqs. (gTSJ) and (g7 

As in the case of the fermion production, we will use the method of Ref. ]I8 to 
compute an approximate expression for the boson occupation number 

n k {r) = 2 sinner), (4.8) 

where the Floquet index or growth factor \i k is determined by: 

cosh(/i fc T) = ReX£ 1) (T), (4.9) 

with x! 1 ' satisfying the same equation (|4.2|) with the initial condition X^(0) = 1, 



L fc 

X^> (0) = 0. The energy density of bosons can be obtained using this approximate 
solution as 



P B (r) = — 2 I dkk 2 n k (r)n k (r). (4.10) 
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Figure 11: The instability chart for bosons coupled to the inflaton (upper panel) and to 
the Higgs (lower panel). Contours correspond to equipotential values of fik = 10~ 3 , 0.01, 
0.03, 0.05, 0.07 and 0.09 in the plane (</ 1(2) , k 2 ). 



and therefore, the ratio of boson to total energies is 

,2 



(2) 



J dkk 2 Vtk{r) sinh 2 (/ifc-r) . 



(4.11) 

Ptotal 7T^i (2 ) ' 

We show in Fig. |TT| the instability chart for bosons coupled to the inflaton (upper 
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Figure 12: The boson growth index /i& as a function of mode number k 2 , for bosons 
coupled to the inflaton. The curves are labelled by the value of its q[ parameter. Note that 
the bands become wider with larger q[. 



panel) and to the Higgs (lower panel). Shaded areas correspond to the instability 
bands for boson production. The darker areas correspond to larger Floquet index 
fXki and thus to a more efficient particle production. In the unshaded areas there is 
no exponential boson production. The instability bands are thinner in the case of 
the coupling to the Higgs compared to those corresponding to the coupling to the 
inflaton, and they shrink to zero for specific values q 2 = 1, 3, 6, 10, etc. This reflects 
the fact that modes with q' 2 = n(n + l)/2 have no instabilities [19| . 

In Fig. [12] we show the Floquet index for bosons coupled to the inflaton and 
different values of the parameters q[ up to 10 6 , corresponding to vertical slices in the 
instability charts of Fig. 11 at those q[ values. The instability bands become wider 
and extend up to larger k 2 with increasing q[ values. Fig. [13] shows the fraction of 
the total energy transferred to bosons coupled to the inflaton. Its mean value grows 
exponentially with time. The exponent is larger for larger q[ values and the boson 
production is more efficient. When this fraction approaches unity, the backreaction 
effect becomes important and the results will be modified ||. These plots have to 
be compared with those in Figs. [| and |H] for fermions coupled to the inflaton and 
Higgs respectively. It is clear that during the first few oscillations of the inflaton, the 
energy transfer to fermions can be more important than that to bosons. 

In the case that the bosons are coupled to the Higgs, as we noticed above, the 
particle production vanishes for q' 2 = n(n + l)/2, and thus no energy is transferred 
to fermions for those q' 2 values. We show in Fig. [H] the fraction of the total energy 
transferred to bosons as a function of time for a set of q' 2 values spanning the range 
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x/T 

Figure 14: The fraction of the total energy transferred to bosons coupled to the Higgs as 
a function of time for different values of the q 2 = 6.5,7.5,8.5,9.5 parameter (labelling the 
curves) . 



between the two zeroes, corresponding to n = 3 and 4. The efficiency of bosonic 
production has an oscillatory behaviour with increasing q 2 . it is maximal for q' 2 ~ 
n 2 /2 and minimal (zero) for q' 2 = n(n + l)/2. 



17 



5. Conclusions 

In this paper we have studied the parametric resonant production of fermions in 
hybrid inflation, with both fields, inflaton and Higgs coherently oscillating after in- 
flation. We have assumed that fermions may couple to either the inflaton or the 
Higgs (or both). The behaviour in the two cases is very different. While fermion 
production is very important in the case of a coupling to the inflaton, even in the 
presence of a bare mass, the production of fermions coupled only to the Higgs is 
generically weak. This is related to the fact that the non-adiabaticity condition, 
dVLk/dr > Ql, is harder to achieve for the Higgs since, when Qk is at a minimum, 
dQ^/dr is also at a minimum (contrary to the inflaton case). When the bare fermion 
mass exceeds the value > y/q~i the fermion production by the coupling to the 
inflaton field is also suppressed. 

We have studied the growth of the fermion energy density and seen that it very 
quickly saturates to an approximately constant value. For fermions coupled to the 
inflaton, the asymptotic value grows with the resonance parameter like h\q\'\ For 
natural values of the couplings, a significant fraction of the inflaton energy can be 
transferred to fermions. On the other hand, for fermions coupled to the Higgs, the 
fermion energy density is of the order p F /ptotai ~ 10 -2 <? 4 /A, which under our working 
assumptions (g 2 <C A) is quite small. 

We have also studied the boson production for both couplings to inflaton and 
Higgs, and compared with the fermion production. While the energy density trans- 
ferred to the parametrically produced fermions saturates after a few oscillations, the 
one in bosons grows exponentially with time. Hence, if both fermions and bosons 
have similar couplings, most of the particle production goes initially into fermions, 
while at late times the boson production is exponentially dominant. There are how- 
ever some values of boson couplings to the Higgs, corresponding to q' 2 = n(n + l)/2, 
which completely inhibit the parametric resonance of bosons. 

When the energy density of the bosonic or fermionic particles produced becomes 
sizeable, they are expected to backreact on the inflaton, affecting its evolution and 
eventually suppressing the parametric production of particles. We have not consid- 
ered this process in detail in this paper since the approach followed is not the most 
suited one for this purpose. For a proper discussion of this issue in the context of 
chaotic inflation, see Refs. |7|. || . 

As a summary, we have shown that the production of fermions in the preheating 
stage of hybrid inflation can be very important. For the range of model parameters 
assumed, hybrid inflation models lead to a more efficient fermion production than 
chaotic inflation models 0, without the need to go to extremely large values of the 
resonance parameter q. Depending on the relative size of the couplings, and on the 
backreaction process, the inflaton energy transferred to fermions may even be larger 
than that transferred to bosons. At any rate, the parametric production of out of 
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equilibrium fermions could have interesting consequences for cosmological issues such 
as the generation of the baryon asymmetry through e.g. the leptogenesis mechanism. 



Note added: while completing the writing of this work a related paper [20 



appeared where the gravitino production during preheating is computed in a super- 
symmetric model of hybrid inflation. 
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